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(—^ , I consider nonleptonic decays of B mesons into two light mesons using the light- 
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cone wave functions for the mesons. In the heavy quark limit, nonfactorizable 
contributions are calculable from first principles in some decay modes. I review 
the idea of the QCD-improved factorization method and discuss the implications 
£^ ' in phenomenology. 

a ' 

\Q ' 1 Introduction 

Nonleptonic decays of B mesons have attracted a lot of attention recently since 
they were observed experimentally in CLEO, BaBar and BELLE. These de- 
cays are important in extracting the information on the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements and CP violation. On the theoretical side, it 
is the least understood area. Since nonleptonic decays involve nonperturbative 
effects such as final-state interactions, it is difficult to obtain a quantitative 
theoretical prediction. However, it has been recently found that nonpertur- 
bative effects such as the nonfactorizable contribution in nonleptonic decays 
could be systematically understood in the heavy quark limit with — > oo. 

Here I will review the current status of understanding on nonleptonic B 
decays very schematically. I will focus on the underlying ideas omitting tech- 
nical complication. 1 hope that this talk will give a clear sketch of what is 
being studied currently. It has been found that nonfactorizable contributions 
in nonleptonic B decays into two light mesons could be calculated using per- 
turbation theory in the heavy quark limit nib — > oo. I will explain in detail 
how we treat nonleptonic B decays in the heavy quark limit and discuss some 
; "*] ■ phenomenological aspects. 

The theoretical framework for B decays is to use the effective weak Hamil- 
tonian at the renormalization scale fi m m&, which is schematically given as 
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where Vij are the CKM matrix elements, Cj are the Wilson coefficients and 
0\ are the effective four-quark operators. The Wilson coefficients are calcu- 
lable order by order in perturbation theory, and the main issue in analyzing 
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nonleptonic decays is how to evaluate the matrix elements of the four-quark 
operators. For example, if B decays into the final two mesons Mi and M2, 
how do we evaluate (MiM 2 \Of\B)? 

First we can use the naive factorization in which we neglect the strong 
interaction effects between mesons and separate the operators into a current- 
current form, and calculate their matrix elements El. Schematically this process 
can be expressed as 

{MxM 2 \Oi\B) » (Afi|J r /t |0)<Af 2 |J' i |B). (2) 

Here is the current operator in the effective Hamiltonian. Each matrix 
element is parameterized by a decay constant or a form factor, which describe 
intrinsic nonperturbative effects. The decay constant and the form factors 
can be obtained from either experiment or other theoretical techniques such 
as QCD sum rules H. 

Unfortunately this naive factorization is unsatisfactory. First of all, there 
is no justification in neglecting the final-state interactions between mesons. 
The argument of color transparency a can be applied in the case of two final 
light mesons, but it should be proved explicitly if the matrix elements can be 
truly factorized including final-state interactions. Secondly, theoretically the 
naive factorization gives an unphysical result. Decay constants and form fac- 
tors are independent of the renormalization scale /i, and the decay amplitude 
has an arbitrary \i dependence through the Wilson coefficients C,*(//). This 
is because we replace the matrix elements, which depend on fj,, by the decay 
constants and the form factors, which are independent of fj,. Therefore the 
resulting amplitude is unphysical. n 

As a remedy to this problem, Ali and Greub □ suggested to calculate 
the radiative corrections of the operators before taking the matrix elements. 
Following this procedure, we can write the matrix element as 

(MiM 2 \Oi\B) «F(M !as )[<Mi|J M |0)<M 2 |J^|S)l + . (3) 

L J tree 

It turns out that the [i dependence in F(/i, a s ) arising from the radiative cor- 
rections cancel the \i dependence in the corresponding Wilson coefficients at 
any given order, hence the decay amplitude does not depend on the renormal- 
ization scale. I call this procedure the "improved factorization" . 

The unphysical \x dependence is absent in the improved factorization, but 
it has other significant problems. First of all, in order to avoid infrared di- 
vergence, F(fj,, a s ) is calculated with off-shell external quarks U. The external 
momenta p 2 play a role of the infrared cutoff. However, because of this, the 
decay amplitudes depend on the choice of the gauge, and the employed reg- 
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ularization schemes a. Therefore the decay amplitude depends on arbitrary 
gauges and the schemes, hence also unphysical. 

If we put external quarks on their mass shell, the gauge dependence and 
the scheme dependence go away, but, in this case, the amplitude becomes 
infrared divergent and gives also an unphysical result. Therefore we need a 
consistent scheme which solves all the problems I mentioned above. Before 

1 explain the main idea of the QCD-improved factorization, note that the 
Feynman diagrams from which infrared divergence appears are those with 
vertex corrections of the weak currents. 

2 QCD-improved factorization 

The source of infrared divergence in the radiative corrections of the effective 
weak Hamiltonian, as mentioned above, suggests an interesting idea. The 
Feynman diagrams which cause infrared divergence are the radiative correc- 
tions of vertices. In other words, the infrared divergence comes from the 
radiative corrections for the decay constant of a light meson and form factors 
for B — > M 2 . 

It reminds us of the hadron-hadron scattering, in which the infrared diver- 
gence of the scattering amplitude is absorbed in the redefinition of the parton 
distribution functions. The remaining hard scattering amplitude and the par- 
ton distribution functions are factorized. We can apply the same idea to B 
decays. The infrared divergences can be attributed to the renormalization 
of the decay constant and the form factors, and other radiative corrections 
constitute nonfactorizable contributions. That is, the infrared divergence is 
absorbed in the definition of the lightjcone meson wave function or the form 
factors. This is first observed in Ref. 

Recently Beneke et al. u formulated this problem in the heavy quark limit 
and extensively studied nonleptonic B decays into two final-state mesons. 
The idea can be summarized as follows: We first take the heavy quark limit 
TO5 — > 00. In this limit we can calculate nonfactorizable contributions system- 
atically in perturbative QCD. Also we can obtain corrections to the pertur- 
bative results by expanding in powers of Aqcd /mi,. 

The next step is to arrange external quarks using Fierz transformation 
such that the quark-antiquark pair, which forms a meson, is included in a sin- 
gle current. And we use the light-cone wave function for the mesons. If we can 
use the operators in the effective Hamiltonian as they are in arranging quarks, 
we call that configuration of quarks as "charge-retention configuration" . If we 
have to switch some quarks using Fierz transformation, we call that config- 
uration as "charge-changing configuration" . This arrangement is important 
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since it determines which radiative corrections correspond to the corrections 
to the decay constants or the form factor. Since the decay amplitude in each 
process depends on how we arrange quarks, the scattering amplitude becomes 
process-dependent. 

The infrared divergence is attributed to the renormalization of the wave 
function or the form factors. And we calculate all the nonfactorizable contri- 
butions along with the spectator contribution and the annihilation channels. 
This procedure is called the " Q CD-improved factorization" . If all the non- 
factorizable contributions are infrared finite, and suppressed as rrib goes to 
infinity, these processes can be treated in the QCD-improved factorization. 
In this case, we have a theoretical method to analyze nonleptonic B decays 
from first principles. In the QCD-improved factorization method, we can 
formally write the matrix element as 



(MiM 2 \O l \B) a (M 1 \J li \0)(M 2 \J"\B) 1 + 0(a s ) + 0(A qcd /to 6 ) 



(4) 



If we neglect the radiative corrections and the 1 / rrib corrections, we restore the 
result obtained in the naive factorization. The matrix element can be written 
as a product of a decay constant and a form factor. And we can systematically 
calculate the corrections to the result in the naive factorization. 

Before we discuss the details of the QCD-improved factorization, I would 
like to explain another approach to study nonleptonic B decays. Keum et al. 
□ have also considered nonleptonic B decays using light-cone wave functions, 
but they concentrated on the calculation of the form factors instead of non- 
factorizable contributions. They calculated a single gluon exchange which is 
responsible for the correction to the form factor. There also appears infrared 
divergence in the form factor, but they introduce the Sudakov factor which 
smears the endpoint region such that there is no infrared divergence. And 
the origin of imaginary parts in their calculation comes from the modification 
of the propagator with transverse momentum, which is totally different from 
the source of the imaginary part in the QCD-improved factorization. 

They do not consider nonfactorizable contributions at the moment, and 
hence the leading-order Wilson coefficients are employed. Because what is 
calculated is totally different in the two approaches, we have to be careful 
when we try to compare the results in both approaches. 

3 Nonfactorizable contribution 

Let us go into the detail of how to calculate nonfactorizable contributions. 
When we use the light-cone meson wave functions for exclusive decays, the 
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transition amplitude of an operator Oi in the weak effective Hamiltonian is 
given by 



{M x M 2 \O t \B) = Vf^ M! f 

Jo 



dxT i I j (x)4> Ml (x) 



+ / d^dxduTi (£,x,u)(I)b(0 < I ) m 1 (x)4>M2( u ), (5) 



where F?^ 2 are the form factors for B — > M 2 , and <pMi(x) is the light-cone 
wave function for the meson Mj. T^(x) and T/ 7 (£,a;,u) are hard-scattering 
amplitudes, which are perturbatively calculable. The second term in Eq. (|J) 
represents spectator contributions. 





Figure 1. Feynman diagrams for nonfactorizable contribution at order a s . The dots repre- 
sents the operators Oi. 



The relevant Feynman diagrams for are shown in Fig. 1. Each Feyn- 
man diagram has an infrared divergence. But if we sum over all the Feynman 
diagrams and symmetrize with respect to x <-> 1 — x, where x is the mo- 
mentum fraction of the outgoing meson, we have an infrared- finite result. 
Another feature of the nonfactorizable contribution is that there appears an 
imaginary part due to the final-state interactions. This plays an important 
role in studying CP violation in nonleptonic decays. The strong phase is cal- 
culable in perturbation theory. And since we are working at next-to-leading 
order accuracy, the dependence on fi becomes very mild. 

There are other nonfactorizable contributions such as the spectator con- 
tributions. The Feynman diagrams for are shown in Fig. 2. In calculating 
nonfactorizable contributions, we use the light-cone wave functions. The pro- 
jection of the quark bilinears to each pseudoscalar light meson wave function 
to the order of twist three can be written as 

(P(p)\q a (y)q' (x)\O) = ^ £ ^ito-v+d-*)*-) 



chay: submitted to World Scientific on February 1, 2008 5 



V V 



Figure 2. Feynman diagrams for the spectator contribution. 



For the B meson, we can write the projection as 



<j>a(u) 
6 




(6) 



(7) 



Here <j> is the leading-twist wave function, and 4> p and ip a are twist-three wave 
functions for the pseudoscalar and the tensor currents respectively. For the B 
wave function, we take the leading-twist wave function only. Since we expand 
in powers of 1/mf,, at leading order 4>b(0 ~ S(£ — Aqcd/to;,). In calculating 
the spectator contribution, we have the integral of the form 



which is enhanced. 

We can consider corrections of order 0{Kqcti /mb). The most important 
contribution comes from the term proportional to jjLp/rrib which is given by 



where mi and ni2 are the masses of the quarks which form a meson with 
mass mp. Compared to other O(Aqcd) terms, this is numerically large. For 
the case of ir + , for instance, it is /j, n + ~ 1.4 GeV. Therefore it has been of 
great interest to calculate higher-twist effects proportional to [ip. However, 
the spectator contribution from 0\ and O2 with higher-twist wave functions 
is infrared divergent 113. It has been suggested that we introduce some pa- 
rameters to regulate the infrared divergence and regard the parameter as a 
theoretical uncertainty. But these contributions are numerically large and 
especially the extraction of the strong phase becomes too ambiguous to say 
anything quantitatively. 
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The annihilation topology poses another problem. When we calculate the 
spectator contribution with the operator O5, the amplitude is also infrared 
divergent. Therefore the annihilation topology can give a significant power 
correction to the decay amplitude. The analysis on the annihilation channel 
with radiative correction is in progress. 

One way to look at the infrared divergence is that it is not a serious 
problem. The divergence comes from the endpoint and it actually gives the 
logarithmic enhancement as 



Since the cutoff Aqcd is actually arbitrary, we can parameterize this con- 
tribution and regard this as a theoretical uncertainty. But its magnitude is 
numerically large, and thus it enlarges theoretical uncertainty. 

Another view is a more conservative one. If there appears an infrared 
divergence in the hard scattering amplitude, the effect of soft gluon exchange 
is really significant and the QCD-improved factorization breaks down at this 
order. It remains to be seen if the QCD-improved factorization really breaks 
down, or there are some other contributions which cancel the infrared diver- 
gence rendering the final result infrared finite. 

4 Conclusion 

The understanding of nonlcptonic B decays into two mesons has acquired a 
new sophisticated level. In the heavy quark limit, nonfactorizable contribu- 
tions are calculable using perturbative QCD for light final-state mesons. When 
one of the final-state meson is heavy, we can still use the QCD-improved fac- 
torization for the case in which the spectator quark in the B meson goes to the 
heavy meson in the final state. If the spectator quark goes to a light meson, 
as in class II decays, the nonfactorizable contribution is infrared divergent, 
and the effect of soft gluon exchange is significant. 

But the analysis of higher-twist effects is yet far from satisfactory. For 
example, the spectator contribution which is proportional to the twist-three 
contribution of the meson wave function is infrared divergent. And the anni- 
hilation topology also has the infrared divergence. The status of the QCD- 
improved factorization method for nonleptonic B decays into two mesons is 
not complete until we disentangle the infrared divergence to give a quantita- 
tive prediction. 

It will be an interesting project to combine the QCD-improved factor- 
ization with the calculation of the form factors using the light-cone wave 
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(10) 



functions. Currently, in the QCD-improved factorization, we use the form 
factors extracted from experiment, as in semileptonic B decays. On the other 
hand, in Ref. El they only consider the calculation of form factors. It will be 
interesting to see whether we can give a consistent theoretical description of 
nonleptonic B decays combining these two approaches. 
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